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We investigate the feasibility of precision frequency metrology with large ion crystals. For clock 
candidates with a negative differential static polarisability, we show that micromotion effects should 
not impede the performance of the clock. Using Lu^ as a specific example, we show that quadrupole 
shifts due to the electric fields from neighbouring ions do not significantly affect clock performance. 
We also show that effects from the tensor polarisability can be effectively managed with a compen¬ 
sation laser at least for a small number of ions (< 10^). These results provide new possibilities for 
ion-based atomic clocks, allowing them to achieve stability levels comparable to neutral atoms in 
optical lattices and a viable path to greater levels of accuracy. 
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I. INTRODUCTION 

The realisation of accurate, stable frequency references 
have enabled important advances in science and technol¬ 
ogy. Well-known examples include the Global Position¬ 
ing System, and tests of fundamental physical theories. 
Increasing levels of accuracy and stability continue to 
be made with atomic clocks based on optical transitions 
in isolated atoms m- By now a number of groups 
have demonstrated superior performance over the cur¬ 
rent caesium frequency standards with the best clocks to 
date having inaccuracy at the 10“^® level mm- Eor the 
past decade single ion clocks have held a leading position. 
However in recent years, advances in laser stability have 
allowed neutral atoms to take advantage of large numbers 
of atoms giving superior performance in stability while 
maintaining some of the best accuracies m- Ion-based 
clocks have been limited to single atoms predominately 
due to the fact that trap influences such as micromotion 
are difficult to control for multiple ions. Indeed, for the 
Al'*' clock, micromotion is a dominant factor in the over¬ 
all error budget [1] . We note that consideration has been 
given to clocks based on small strings of less than 10 ions 
stored in linear radio-frequency (RE) traps |11| . 

Micromotion is caused by the RE drive used to con¬ 
fine the ion. It causes two correlated effects: the rapid 
oscillatory motion at the RE frequency H gives a second- 
order Doppler shift, and the electric field driving the mo¬ 
tion induces an AC stark shift proportional to Aa where 
Aq: = ae — ctg is the differential static scalar polarisabil¬ 
ity. In 1998 [H] it was noted that, when Aa < 0, these 
two effects could be made to exactly cancel for a well 
chosen value of D. This value, which we refer to here as 
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the magic RF frequency in analogy with the magic wave- 
length for optical lattices, depends only on the properties 
of the atom. The existence of such a magic RF frequency 
gives rise to the important question of what happens in 
a large ion crystal where micromotion effects can be very 
pronounced at the edges of the crystal. If the cancelation 
is maintained, then precision frequency metrology with 
large ion crystals is potentially feasible. 

In this paper we show that dominant higher order mi¬ 
cromotion effects can also be mitigated, and that clock 
operation would not be affected by micromotion. For 
atoms with an electronic angular momentum J > 1/2, 
we show that quadrupole shifts due to the electric fields 
from neighbouring ions do not significantly affect clock 
performance. We also show that shifts arising from the 
tensor polarisability can be effectively compensated with 
an additional laser field. Together these results show that 
ions having Aa < 0 can reap the benefits of large num¬ 
bers of ions, just as neutral atoms can in optical lattices. 
We illustrate our analysis using ^’^®Lu+ as a concrete ex¬ 
ample m. but the ideas can be readily adapted to any 
other ion with Aa < 0. 


II. MICROMOTION 

We start by writing the RF and static electric field 
potentials in the form 

, mUlujz T. ^ , TOwJ rp 

</rf = -r ArfrcosDt, A^r (1) 

2q 2q 

where is one of the pseudo-potential oscillation fre¬ 
quencies, Ul is the RF drive frequency, and m and q are 
the mass and charge of the ion respectively. The matrices 
Arf and As determine the curvatures of the potentials and 
in general, we may choose A^f to be diagonal. Defined in 
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this way, the pseudo-potential approximation is 

V (r) = r (2) 

If we scale time by 2/Vl and length by 

J ’ 

then the equations of motion (e.o.m.) are given by 

Vi + (e^ Ag -I- 2eAi.f cos 2t) Vi — ^ ^ = 0 (3) 

where e = 2u3^l^. Following the treatment in [2] we 
assume a stable 7r-periodic crystal solution exists, which 
may be expressed as a Fourier expansion 


at which micro-motion shifts cancel as first pointed out in 
|12) . For one ion, Eq.j^is sufficient for even the very best 
clocks [DEIS]. However, for large ion crystals, higher 
order corrections should be considered. To this purpose, 
we first note that the n = 2 Fourier component of the 
e.o.m. is 

R-4,i = T—ArfR2,i, (H) 
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to lowest order. Substitution into Eq. j^then gives 
^Ag -f — A^f^ — 41^ R2,i + eAi-fRo.i 

-e"I]Qb(R2,*-R2j) =0. (12) 

Using the fact that (l — ^ « I -|- e^A, we can solve 

for R 2 y to get 


vUt)= Y. (4) 

n— — oo 


In this form, Rgy is the time-averaged position of the 
ion. Substituting this expansion into the e.o.m. and 
using a Taylor series expansion for the Coulomb term 
about Ro.ij = Ro,i — Ro,i we can obtain an infinite set 
of coupled equations for R 2 ri,i representing the Fourier 
expansion of the e.o.m. Defining 
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the first two Fourier equations are given by 


e^AgRo i + 2eArfR2,i — 'Y, ^ (®) 


and 


(e^Ag — 4l) R2_i + eArf (Ro,i -f R4,i) 

-e"I]Qb(R2.*-R2.j) =0. (7) 

To lowest order, Eq.[^ gives 

R2,z = ^ArfRo,i (8) 

which expresses the fact that the micro-motion ampli¬ 
tude is directly proportional to the RF electric field at 
the position of the ion. With this approximation, the 
fractional shift of a clock transition is given by 
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( 9 ) 

For Aa < 0 this leads to a magic RF drive frequency 
defined by 
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( 10 ) 


R2,. = ^ (^Ag + eArfRo,, 

— Yq Qij^rl (Ro,i — Roj ) • (13) 

3¥=i 

The term on the second line of Eq. is the coupling 
of the micromotion amplitudes of each ion through the 
Coulomb interaction. Physically it arises from a distor¬ 
tion of the space-charge potential due to differential mi¬ 
cromotion amplitudes between ions, which provides an 
effective RF electric field in addition to the trap drive. 
This effective electric field also provides an additional AC 
stark shift. 

The electric field at the ion due to the space charge 
is 

E, = ^ (F„- + 2Q,, (R 2 .* - R 2 ,,) cos(2t)). (14) 

^ 3¥=i 

Using Eq. I^and defining 

Woy = Y^ QbArf(Roy — Roj); (15) 


the amplitude of the net RF electric field on the ion 
is then 


Ej,rf = — - 
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Eq. 16 and noting that (l-fe^A)" 


Using Eq. 

I-|-2e^A the fractional frequency shift of a clock transition 
due to terms oscillating at the RF drive frequency is then 


Avj 

V 


ujj 




RqiArfRo.i 


, n 

% 


Ro_iArfWo,i 
1 


“k 2 ( Ag + 20 '^i'f ) ArfRo,i f • (17) 


















3 


To the same order of approximation, we must include a 
DC stark shift from the space charge and a time dilation 
shift from the higher harmonic R 4 ,i. These can be added 
independently, which follows from the orthogonality of 
the Fourier components. At equilibrium, the sum of all 
DC fields on the z*** ion exactly balances the pseudo¬ 
potential force from the RF field so we have 


Ej,dc = V 


1 mujJ ,2 


ArfR-o,i- 


2 q 

This gives a fractional frequency shift 


(18) 


Avi Aa / muj'^V 




and the time dilation shift from R 4 i is 


Aiy, _ eV ujJ\ T « 4 t } 

RoAfRo,.. 


(19) 


( 20 ) 


The shifts given in Eq. [T^ & can then be added to 
Eq. 
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to give the total shift correct to 0{e 


In addition to a fractional frequency shift, micromotion 
also gives rise to an effective frequency modulation of the 
probe; a fact that can be used to detect micromotion m- 
The effect reduces the probe coupling to Jo(/3) where Jq is 
the first order Bessel function and the modulation index 
is given by /3 = 21k • R 2 ,i. Because the micromotion 
amplitude can be very large for ions far removed from 
the zero point of the RF field [T^], the ability of the 
probe to drive the clock transition can be significantly 
diminished. This can be avoided by probing along the 
RF null axis of a linear Paul trap, for which there is very 
little micromotion. 

For the linear Paul trap, we can write 


Arf = 
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where a determines the strength of the RF conhnement 
relative to the DC held. The parameter 6 determines 
the asymmetry in the transverse dimension and can be 
tuned arbitrarily by appropriate choice of biasing volt¬ 
ages. Writing A^f = oA, the total shift from Eq. [T^ [T^ 
& [ 20 ] can be written 
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where the Afc are given by 
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The hrst term of Eq. 22 is zero for fl = Dq-Ao and this 


can be viewed as a correction to the magic RF frequency 
given by Eq. It also applies to the single ion case 
even though this term includes the shift due to the DC 
component of the space charge. This is because it is 
the space charge that provides the static electric held 
necessary to displace an ion to the equilibrium position 
Roy. The second term also applies to the single ion case 
and is only present when the transverse conhnement is 
non-degenerate. The third term is only applicable to the 
many ion case as it is a consequence of the induced RF 
held from the oscillating space charge. 

When probing a clock transition, micromotion will give 
rise to an inhomogeneous broadening of the line and a 
shift of the line centre. Since the hrst term in Eq. 
can be tuned to zero, the inhomogeneous broadening is 
limited only by the remaining terms and the degree of 
broadening is determined by the size of the crystal. The 
second term is suppressed for small S and, at 11 = floVAj 
the hnal term scales as . If we can neglect the inhomo¬ 
geneous broadening, the shift of the line centre is then 
determined by the average of Eq. For large crys¬ 
tals we can assume the density of ions is approximately 
constant and the average of Eq. can be taken as an 
integral over a continuum. The shift of the clock tran¬ 
sition then scales as where N is the number ions 

and the scale factor depends on the geometry of the trap 
|I5) . This leads to a further modification to the magic 
RF drive frequency at which the scale factor vanishes. In 
a practical application we would simply vary the number 
of ions and tune the variation of the clock frequency to 
zero by adjustment of D. 

For reasons to be discussed in the next section, it is 
advantageous to take a spherically symmetric trap for 
which a = 1/3 and 5 = 0. In the spherically symmet¬ 
ric case, an analytic approximation for Woy can be ob¬ 
tained. By taking the continuum limit with a constant 
density of ions, we can approximate the sum in Eq. 
by an integral and we obtain 
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(24) 


In this case the fractional frequency shift takes the simple 
form 
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and a; = 1 + 


(26) 


where the A), are given by 
A'-l 

Note that Aq — Ag ~ 0{e^) indicating that the oscillating 
space charge has no significant affect to the accuracy of 
the treatment given. 

At this point it is useful to illustrate our analysis with 
an example. We simulate the distribution of ions by in¬ 
tegrating the e.o.m. in the pseudo-potential approxima¬ 
tion including a small damping term to anneal the initial 
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FIG. 1. (a) Level structure of ^^^Lu"*", which has a nuclear 
spin 1 = 7 , showing the clock transition at 848 nm and the 
cooling and detection transition at 646 nm. A laser at 615 nm 
provides a magic wavelength to compensate tensor polaris- 
ability shifts from the RF field, (b) Transitions need for hy- 
perfine averaging, (c) Transitions needed to effectively realise 
the transitions in (b) due to constraints on the quantisation 
axis - see Sec. mm 


state. For the initial state we use multi-layered Mackay 
icosahedra m with a random offset to each particle co¬ 
ordinate of about 10% of the initial minimum particle 
spacing. With N = 100 to 5000 ions we were able to 
confirm the results given in m- In particular, to a good 
approximation we have 

(Rj^jA^Ro,*) = - 0.3964. (27) 

For the ion properties we use which has a level 

structure as illustrated in Fig. [I] For ^^®Lu+, the esti¬ 
mated value of the differential static scalar polarisability 
is Aa = —2.19 a.u. giving a magic RF frequency 

of flo ~ 27r X 23.2 MHz [TH]. We take Wz = 27r x 200 kHz 
which, for ^’^®Lu“'", gives I = 7.94 fim. So, for the param¬ 
eters given, we have 

^ « 8.3 X 10-1®, « 3.0 x lO""!. (28) 

With H = Hq and N = 5000 ions, the higher order terms 
result in a broadening of 3.3 x and an average shift 
of 1.4 X IQ-i'’. These values would scale as When 

H = Ho-y/A)) the average shift vanishes and this is only 
a 7 X 10“® fractional change to the RF drive frequency. 
Comparison of the clock frequency at different numbers 
would therefore provide an accurate assessment of H and 
Aa. 

From the above analysis it would appear micromotion 
would not limit the accuracy that can be achieved. More¬ 
over, the small broadening effects would not affect achiev¬ 
able stability in the foreseeable future. As the crystal size 
increases, so does the demands on the pointing stability 



of the probe but, for 10^ ions and an angular misalign¬ 
ment of 0.1°, the Rabi frequency for the outermost ions 
would be diminished by just 3 x 10“^. 

Our analysis has only considered the scalar polarisabil¬ 
ity. Levels with J = 1/2 have a vector polarisability but 
this is only relevant for circular polarisations. For candi¬ 
dates such as Lu'*', which have clock states with J > 1/2, 
we must also consider effects arising from the quadrupole 
moment and the tensor polarisability. 


III. CONSIDERATIONS FOR J > 1/2 

In [13] , Barrett showed that averaging over transitions 
to all hyperfine states of a fixed mp < I — J cancels dom¬ 
inant magnetic field effects and quadrupole shifts when¬ 
ever the nuclear spin, I, is at least as large as J. This 
averaging, which we shall refer to as hyperfine averaging, 
is very general and applies to any perturbation that can 
be described by a rank k > 0 tensor operator that does 
not depend on I. Hence, we need only consider the in¬ 
homogeneous broadening arising from such interactions. 
The following considerations can also be readily adapted 
for candidates with 1 = 0, such as ®®Sr+ for example |1|, 
which utilise averaging over all mj. 


A. Quadrupole Shifts 


Since the quadrupole shift from the DC trapping field 
is fixed, it does not contribute to any broadening and so 
we only consider the quadrupole fields arising from the 
space charge. Since the quadrupole field from neighbour¬ 
ing ions falls off cubicly with distance, the quadrupole 
field experienced by an ion is due mostly to its local en¬ 
vironment. Except for ions near to the edge of the crystal 
we can anticipate that the local environment is essentially 
the same for each ion and the resulting quadruple shift 
is relatively homogeneous. Indeed for sufficiently large 
numbers, it is known that a crystal of long range order 
forms [ISlEniEI] and the limiting structure is that of 
a body-centered cubic (bcc) lattice. In this regime, the 
quadrupole shift would be constant for the bulk of the 
crystal and by symmetry it would be zero. 

The tensor describing the quadrupole field for the 
ion due to all other ions is given by 


VEf) = 


e 


3¥=i 


(29) 


where Qij is given in Eq. Itano [52] has derived the 
quadrupole shift for a general orientation of the quanti¬ 
sation axis relative to the coordinate system for which Qi 
is given. This shift factors into a geometric term, a state 
dependent scale factor Cp^mp on the order of unity, and 
an overall scale factor quantifying the size of the shift. 
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The geometrical factor is given by 

(3 cos^ /3 - l) + i sin(2/?) (Qa:z cos a + Qy^ sin a) 
+ I3{{Q XX - Qyy) cos{2a) +2QxySin{2a)), (30) 

where we have used the Euler angle definitions in [22| and 
dropped the subscript i for convenience. For the mp = 0 
states of the ^Di level of ^^®Lu~'', the state dependent 
scale factors are —2/5,1, —3/5 for the F = 6 , 7,8 hyper- 
fine levels, respectively. In the calculations that follow we 
omit this factor. The overall scale factor depends on the 
quadrupole moment. For the estimated value of — 1.3eaQ 
m its size is —1.3mw^aQ « /i x 2.5 Hz where Oq is the 
Bohr radius. 

In Fig. we plot the distribution of quadrupole shifts 
for N = 5000 for a spherically symmetric trap. As in 
we do see a dependence of the final crystal config¬ 
uration on the initial condition for larger JV. Configu¬ 
rations starting from a bcc-lattice tend to stay in this 
configuration with some rounding near the boundary of 
the crystal. Hence we give the distribution of quadrupole 
shifts for two types of initial conditions: a multi-layered 
Mackay icosahedra applicable to smaller numbers of ions, 
and a bcc lattice applicable in the limit of large N. The 
distribution on the left is more applicable to smaller num¬ 
bers as considered here and is approximately Gaussian 
with a standard deviation of 0.078 Hz. The distribution 
does not depend on JV and we have verified that there 
is no significant dependence on the field orientation as 
expected from the spherical symmetry. This level of the 
broadening should not be an issue even for the very best 
lasers available today. 

The distribution on the right is applicable to larger 
numbers of ions. However crystals of long range order 
have been reported for smaller numbers [20) . Compar¬ 
ing the two distributions we clearly see the effect of the 
bcc-lattice component giving the expected peaking of the 
distribution near zero. This distribution does have a de¬ 
pendence on the orientation of the B-field and we have 
oriented the field along one of the axes of the interior 
bcc structure which coincides with the trap axis [53] . We 
would expect the dependence on the orientation of the B- 
field to diminish for larger numbers as the bcc component 
becomes much more prominent. 

The distribution of quadrupole shifts depends only 
on the geometry of the crystal and not on its overall 
size, at least for the range of numbers we have explored. 
Prolate ellipsoidal crystals in more conventional linear 
Paul traps, in which the transverse confinement is much 
stronger than the axial confinement, have a much broader 
distribution of quadrupole shifts. Moreover, the width of 
the distribution depends on the orientation of the trap 
relative to the quantisation axis as may be expected. For 
these reasons, we have restricted our attention to a spher¬ 
ical geometry. 
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FIG. 2. Distribution of quadrupole shifts for 5000 ions in a 
spherically symmetric trap. Initial starting distributions are 
a multi-layered Mackay Icosahedra (left) and body-centered 
cubic lattice (right). 


B. Tensor Polarisability 


The tensor polarisability also gives rise to a shift of the 
clock frequency from the RF fields. As shown in |22l [24] , 
this contribution is given by 


Sp 

V 


Cp^mp Ct2,J 
4 hv 


{iEl 


E^), 


(31) 


where Cp^mp is a state dependent scale factor identical 
to those for the quadrupole shift, (•) indicates a time 
average over one cycle of the oscillating field, and the 
tensor polarisability 02 ,j is in general frequency depen¬ 
dent. Since the RF frequency is small relative to any 
optical frequency of interest, we can use the DC value 
of 02 ,j ~ —5.0 a.u. for the tensor polarisability m- 
If the quantisation axis is aligned along the trap axis, 
then Ez = 0 and the shift has the same form as dis¬ 
cussed for the scalar polarisability. In general, due to the 
F dependent scale factors, we cannot simply modify the 
scalar polarisability to account for the effect [25]. Since 
Eq. [^ applies at all frequencies, we can use a laser field 
to reduce the broadening that arises. We also note that 
the constraint on the quantisation axis forces us to re¬ 
alise mp = 0 to mp = 0 transitions as an average over 
mp = ±1 to mp = 0 as illustrated in Fig. This in¬ 
troduces further averaging to that discussed in [T3| but 
does not affect any of the points considered here. 

Use of a laser to reduce the broadening arising from 
the tensor polarisability requires the spatial dependence 
of the beam to match the spatial variation of the RF field 
and the use of a magic wavelength at which the dynamic 
differential scalar polarisability of the clock transition is 
zero. For Lu'*', such a wavelength can be found for a laser 
tuned between the ^Di to ^Pq and ^Di to ^Pi transitions. 
From matrix elements given in |26j , we find a magic wave¬ 
length at « 615 nm with 02 , j ~ 100 a.u. The sign of a 2 j 
relative to the DC value requires the use of a doughnut 
mode m which, to lowest order, has an intensity profile 
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Uncompensated Shifts 



Compensated Shifts 
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FIG. 3. Distribution of shifts due to the tensor polarisabil- 
ity for N = 1000 ions. Distribution on the right has been 
compensated using a doughnut Laguerre-Gauss beam with a 
waist of 100/ ~ 800 /rm. We have used the same spherical 
geometry as in the previous sections. 


that has the same quadratic dependence on the distance 
from the trap axis as the amplitude of the RF field. 
In this case, off-resonant scattering from the compensa¬ 
tion beam is completely determined by the amount of 
compensation needed. For a crystal of 10"^ ions, we esti¬ 
mate a scattering rate « 0.005 s“^ for the outermost ions. 
Hence, for Lu'*', off-resonant scattering will not limit this 
approach. A more practical limit would be the mode 
matching of the laser profile to the spatial variation of 
the RF field. But we emphasise that the mode match¬ 
ing need only be sufficient to reduce the broadening since 
hyperfine averaging eliminates any residual shift. 

In Fig. we illustrate the compensation of broaden¬ 
ing for N = 1000 ions using a doughnut Laguerre-Gauss 
beam with a waist of lOOZ « 800 fim with the same spher¬ 
ical trap geometry used in the previous sections. Due 
to the fact that the width of the broadening scales as 
it would likely become impractical to go beyond 
N = 1000 ions. This assessment is based on due con¬ 
sideration of demonstrated mode purity of higher order 
Laguerre-Gauss beams IW, and realistic constraints on 
beam size. 


IV. PROSPECTS FOR LU+ 

As a clock candidate, Lu’*' has a number of favourable 
properties leading to low systematic shifts which are sum¬ 
marised in Table |lj The blackbody radiation shift at 
300 K is based on the current estimate of the differential 
static polarisability of Aa = —2.19 a.u. HZ]. The second- 
order Doppler shift due to residual secular thermal mo¬ 
tion assumes Doppler limited cooling on the 646 nm tran¬ 
sition. This has a line-width of F = 27r x 2.45 MHz pro¬ 
viding one of the lowest Doppler cooling limits amongst 
the ions and yet sufficiently large to allow a collection 
of > 5 photons/ms per ion during detection. Hyperfine 
averaging cancels dominate Zeeman shifts leaving only a 
residual quadratic shift of « 5Hz/mT^ due to coupling 


Effect 

Shift (10-1®) 

Blackbody Radiation @ 300K 

53.3 

Secular Doppler 

-0.05 

Micromotion 

0 

Quadrupole Shifts 

0 

Quadratic Zeeman @10 /rT 

-1.4 

Probe AG Stark (200 ms vr-pulse) 

< 50 


TABLE I. Summary of systematic shifts for Lu’^ after hyper¬ 
fine averaging. 


to the ^D 2 fine structure level [T3j. The shift given in 
the table assumes an operating field of 10 ^T. Finally, 
AC Stark shifts from the probe laser are based on cur¬ 
rent estimates of the dynamic polarisabilities at the clock 
frequency (Aao = —18.12 a.u. and a 2 = —12.44 a.u.), 
and a lifetime of 62 h HZ]. The inequality given in the ta¬ 
ble applies to all six transitions involved in the hyperfine 
averaging as indicated in Fig. [^c). 

We can expect much of the AC Stark shift from 
the clock laser to be eliminated by hyper-Ramsey spec¬ 
troscopy |28l [29]. Thus a temperature inaccuracy of 1 
degree at room temperature would permit fractional in¬ 
accuracies below 10“^®. Within a cryogenic environment 
we can also anticipate inaccuracies beyond 10~^®. 

As discussed in Sec. jlllB] we can expect to be lim¬ 
ited in practice to A « 1000 ions. For this many ions, if 
we combine the quadrupole shifts and the compensated 
shifts from the tensor polarisability, we obtain a reason¬ 
ably symmetric distribution which can be roughly ap¬ 
proximated to a Gaussian with a standard deviation of 
about 0.1 Hz. With the actual distribution, a simulated 
Ramsey experiment with a Tm = Is free precession time 
yields an 80% contrast in the Ramsey fringes. For Lu+ 
state preparation and detection can be expected to take 
~ 1% of the total interrogation time and hence the Dick 
effect [30] should not have a significant role. Neglect¬ 
ing the slight loss in fringe contrast, gives an estimated 
projection noise limited stability [51 [5T] of 

/ N 1 1-5 X 10”^^ 

^- /.rr^ - -7=-■ (32) 

2TTVo^NTmT yjT 

With a stability given by Eq. |32[ measurement of 
the clock frequency at levels of 10“^could be achieved 
within K, 5 minutes. Assessment of micromotion shifts 
associated with small inaccuracies of the magic RF fre¬ 
quency does, however, require comparison of clock mea¬ 
surements with different numbers of ions. Integration 
times of approximately 1 day would permit assessment 
of the clock at the 2 x 10“^® level for 100 ions. Sub¬ 
sequent comparison with 1000 ions would then provide 
a measurement accuracy of Aa and the magic RF fre¬ 
quency at the 10“® level. Using micromotion shifts to 
determine Aa has been demonstrated with a single ion 
[32] . With many ions, the sensitivity of this approach is 
substantially improved. 
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V. ALTERNATIVE TRAPS 


Our approach has focused on linear Paul traps. Other 
approaches may be feasible such as multipole traps or 
Penning traps. In multipole traps low numbers of ions 
initially populate a single 1 -dimensional ring of ions 
that can hold up to several tens of ions. By symmetry 
quadrupole shifts and shifts from the tensor polarisabil- 
ity are practically constant for all ions. As more ions 
are added, more rings form. For more than 2 rings, the 
shifts split into multiple values that can be separated at 
the Hertz level. Thus it would be difficult to go beyond 
a few 10^ ions by this approach. Nevertheless, this could 
be achieved with very little broadening thus allowing for 
much longer Ramsey times. 

Penning traps have been used to confine and control 
very large numbers of ions |21j . Ion confinement in a 
Penning trap is due to the ion crystal rotation through 
a large uniform magnetic field. In a Penning trap the 
combined fractional frequency shift is 

1 / 2Vc/ hiy \ e ) 

where p is the cylindrical radius of an ion in the crys¬ 
tal and ujr is the rotation frequency of the crystal. This 
then leads to a magic rotation frequency analogous to 
the magic RF frequency for the RF Paul trap. For the 
Penning trap there are no higher-order corrections so 
the 2nd order Doppler and polarization compensation 
should work very well, but there are constraints on uJr 
due to available magnetic fields. Specifically, the rota¬ 
tion frequency is bounded by the cyclotron frequency 
Dc = eBjm |33]. This leads to the constraint 



B > 




hv 

—c^Aa 


(34) 


Unfortunately this favours a large polarisability and 
small clock frequency. Even a clock frequency as low 
as 10^^ Hz and a differential polarisability Aa = —100 
requires B > 22 T. 


VI. CONCLUSION 

We have shown that ions with a negative differential 
static polarisability should allow high precision metrol¬ 
ogy on large ion crystals. More specifically, we have 
shown that micromotion does not give rise to any signif¬ 
icant inhomogeneous broadening and that higher order 
frequency shifts can be managed through adjustment of 
the magic RF frequency. For clock candidates that sup¬ 
port a quadrupole moment, we have shown that spher¬ 
ically symmetric traps show very little broadening due 
to quadrupole shifts induced by neighbouring ions and 
this broadening is not dependent on the size of the crys¬ 
tal. We have also shown that broadening arising from 
the tensor polarisability can be compensated by a laser 


field, at least for smaller numbers of ions (< 1000). This 
will extend the advantage of using large numbers to ion- 
candidates having Aa < 0; an advantage that has al¬ 
lowed neutral atoms to surpass the performance of single 
ion standards. In the case of ^^®Lu“'', this approach could 
outperform the current state of the art by an order of 
magnitude in both stability and accuracy. 

We have not included effects due to anharmonicities of 
the trapping fields as these effects are design dependent. 
However the framework we have used follows that given 
in m and should allow such effects to be included given 
a particular design. The effect will be to give a spatial 
dependence to Ag and A^f. But the influence of Ag only 
appears to second order in e and we have shown that 
these effects do not contribute significantly under most 
circumstances. Furthermore, the spatial dependence of 
Arf in Eq. would not change the lowest-order equation 
for the magic RE frequency. Hence, we believe that the 
main effect of anharmonicity will be to affect the higher 
order terms only. This would introduce a small amount 
of broadening and not change the general conclusions we 
have made here. We may anticipate a variation in crys¬ 
tal density giving further broadening due to quadrupole 
shifts, but this would likely only be significant for highly 
anharmonic confinement. 

We have also not considered magnetic field inhomo¬ 
geneities as these are again design dependent. Magnetic 
fields arising from currents induced by the trap driving 
field [T] can be expected to have a significant spatial vari¬ 
ation giving rise to a broadening through the quadratic 
Zeeman shift of the clock states. Static field inhomo¬ 
geneities would also give rise to additional broadening 
through the linear Zeeman effect. Thus candidates with 
small B-field sensitivities would be desirable, but this 
is true of any clock. Moreover, it may be possible in 
a specific implementation to compensate any significant 
broadening with additional fields as we have shown for 
tensor polarisability effects. 

Although the importance of Aa < 0 was pointed out 
over fifteen years ago [ 12 ], it has not played a signifi¬ 
cant role in the development of ion-based atomic clocks. 
This is perhaps due to the scarcity of candidates having 
this property. To our knowledge there are eight can¬ 
didates that have been reported in the literature: B+ 
[all , Ca+ [all, Sr+ [an, Ba+ [35], Ra+ [35], Er2+[T9], 
Tm^+[T2], and Lu'''[T2]. Of these candidates, B+ is the 
only candidate with aJ = 0toJ = 0 clock transition 
for which quadrupole and tensor polarisability restric¬ 
tions do not apply. However, the magic RE frequency for 
B+ is « 800 MHz which may be technically challenging 
to implement. Moreover, the only cooling and detection 
channel available is the ^Sq to ^Pi at 137nm. 

From the alkaline-earth metals, Ba'*' is an interesting 
possibility. For ^^^Ba’*', there are a number of states 
of the D 3/2 level with Cp^mF — 0- Hence, quadrupole 
and tensor polarisability considerations would not apply. 
Since these are the main limitations to working with large 
numbers, very high levels of stability could be possible 







with this ion. It is hoped that this discussion spurs in- ACKNOWLEDGMENTS 
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